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ABSTRACT 

We present uniform CFHT Megacam g and r photometry for 34 X-ray selected galaxy 
clusters drawn from the X-ray Multi-Mirror (XMM) Large Scale Structure (LSS) sur- 
vey and the Canadian Cluster Comparison Project (CCCP). The clusters possess well 
determined X-ray temperatures spanning the range 1 < fcT(keV) < 12. In addition, 
the clusters occupy a relatively narrow redshift interval (0.15 < z < 0.41) in order 
to minimize any redshift dependent photometric effects. We investigate the colour bi- 
modality of the cluster galaxy populations and compute blue fractions using criteria 
derived from Butcher and Oemler (1984). We identify a trend to observe increasing 
blue fraction versus redshift in common with numerous previous studies of cluster 
galaxy populations. However, in addition we identify an environmental dependence of 
cluster blue fraction in that cool (low mass) clusters display higher blue fractions than 
hotter (higher mass) clusters. Finally, we tentatively identify a small excess population 
of extremely blue galaxies in the coolest X-ray clusters (essentially massive groups) 
and note that these may be the signature of actively star bursting galaxies driven by 
galaxy- galaxy interactions in the group environment. 



, 1 GALAXY POPULATIONS IN CLUSTERS 

' The observed properties of galaxy populations reflect the en- 
vironment in which they are located. Comparisons of galaxy 
, populations drawn from low (the fleld) and high density 
' (rich galajcy clusters) environments indicate that the popu- 
lation distribution descr ibed using me asures such a s current 
star f ormation rate (e.g. iBalogh et al] 1999, Poggia nti et al.l 
l2006 h. integr ated colour (e.g . Blanton et al. 2005) and mor- 
phology (e.g. iDressler et al.l il997. Trcu et al. 2003)) varies as 
a function of changing environment. From studies such as 
these it is clear that galaxies located in the cores of rich clus- 
ters display lower star formation rates, redder colours and 
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more bulge dominated morphologies compared to galaxies 
located in the field. 

Studies of the fraction of blue galaxies contributing to 
a galaxy cluster provided some of the first direct evidence 
for the p hysical transformation of galaxies in cluster envi- 
ronments. iButcher fc Oemleil l|l984h reported an increase in 
the fraction of blue galaxies in 33 rich galaxy clusters out 
to 2 ~ 0.5 compared to local clusters. However, subsequent 
studies designed to expand upon this initial discovery high- 
lighted the many complexities associated with this relatively 
straightforward technique, including varying intrinsic cluster 
pr operties with redshift (e. g. X-ray luminosity as discussed 
bv lAndreon fc Ettori|[l999l ). the use of fc-corroctions to de- 
termine rest frame colour distributions (lAndreon 200 51 and 
the challenge of obtaining large sa mples with unif orm pho- 
tometry (successfully overcome bv lLoh et allbOOSl ). 

In addition to the Butcher-Oemler (BO) effect mea- 
sured employing optical galaxy colours, analogous BO-type 
effects have been reported as either a morphological BO ef- 
fect (inc reasing spiral fraction in clusters with increasing 
redshift; IPoggianti et al.l Il999l ) and an infra-red (IR) BO 
effect (increasing fraction of dust enshrou ded star forming 
falaxies in clusters wi th increasing redshift: IDuc et aL l l2ooa 
ISaintonge et al.ll2008l ) to name two examples. 
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An alternative approach is to consider an environmen- 
tal Butcher-Oemler effect whereby one attempts to deter- 
mine the variation of blue fraction as a function of vary- 
ing intrinsic cluster properties selected over narrow redshift 
intervals. This has been achieved by comparing blue frac- 
tions within clusters at increasing c lustercentric radii (e.g. 
lEllingson et al.ll200ll : lLdr et al I2OO8I ) or by considering blue 
fractions me asured bet ween clusters of differing X-ray lumi- 
nosities fe.g. IWake et al. 2005). A number of studies report 
the decrease of the fraction of blue galaxies with decreas- 
ing scaled clustercentric radius, e.g. the viri al radius deter- 
mine d employing either clust er dynamics (|Ellingson et all 
correlation properties (ILoh et al. 1|20QJ), or extrapo- 
lated from X-ray properties (|Wake et al.l 2005f ) . 

The currently favoured explanation for these observed 
trends is that infalling field galaxies are processed physi- 
cally as they travel from the field, through th e cluster out- 
skirts and virialise in the central cluster region l|Berrier et al.l 
I2OO9I ). However, the extent to which cluster galaxies were 
"pre-processed" by physical effects occuring durin g an ear- 
lier r esidence in a galaxy grouj^] remains contested (|Li et all 
I2OO9I ). Numerous physical processes have been suggested as 
the agents of this apparent transformation of galaxy popula- 
tions. However, the dominant physical process(es) to which 
an infalling galaxy is subject remains unclear. 

The two principal examples of such processes are ram 
pressure stripping and galaxy galaxy interactions. Ram 
pressure stripping describes the effective force experienced 
by the diffuse gas component of the infalling galaxy as 
it tra vels through the ho t, dense intra-cluster medium 
(ICM; iGunn fc GottI 11972*). Both hot and cold gas may 
be stripped from the infalling galaxy leading to the ex- 
haustion of the av ailable gas supply that will otherwise 
cool a nd form stars (|Abadi et al.]|l999l : iKawata fc MulchaevI 
l2008l : iMcCarthv et all I2OO8I '). Hvdrodvnamical simulations 
have indicated that the ICM associated with galaxy group 
and cluster environments will strip of order 70% to 100% 
respectively of hot gas from a typical infalling spiral 
galaxy within on e crossing time iKawata fc Mulchaevlliooi : 
[McCart hv et al.l 12008). The effect is manifest as a sharp 
decline in the galaxy star formation rate effective on a 
timescale comparable to the rate at which the unstripped 
cold gas supply is consumed. The observations of a popula- 
tion of red (i.e. passive) spiral galaxies contributing to the 
red sequence would appear to support the view that some 
fraction of galaxies e xperience the stri pping of disc gas via 
ram pressure effects (|Wolf et all [20091 ). The observation of 
:alaxies displaying extended Hi tails in the Virgo cluster 
Chung et al.ll2007l ) is nominally consistent with the expec- 
tations of ram pressure stripping. However, we note that the 
authors also report the presence of close companions to these 
galaxies in a number of cases and comment that a combi- 
nation of ram pressure and tidal stripping provides a more 
compelling explanation. 

When referring to galaxy-galaxy interactions we note 
that this may indicate one of a wide range of encounters. In- 
teractions may be predominantly tidal between close neigh- 



^ Note that we generally refer both groups and clusters of galaxies 
as "clusters" in this paper and quantify our description using the 
X-ray temperature. 



bours resulting in halo gas being moved outwa rd where it is 
more readily stripped (e.g. IChung et al.ll2007l '). High speed 
encounters (either referred to as harrassment or threshing in 
the literature) may also result in tidal stripping of halo gas 
jMoore et al.|[l99a ). Finally, infalling galaxies may merge 
with existing cluster members. The products of such merger 
encounters may be predicted by considering the mass ratio 
of the merging galaxies: large mass ratios result in enhanced 
star formation in the satellite galaxy yet may not lead to a 
star burst in the more massive companion l|Cox et al.ll2008l ). 
Equal mass mergers on the other hand result in the com- 
plete disruption of the infalling spiral galaxy (for example) 
to form a bulge dominate d system ac companied by a sig- 
nificant central star burst l|Di Matteo et al. ,200'^ ). The in- 
ternal disruption associated with such strong interactions 
can result in a short term enhancement in star formation 
followed by a rapid exhaustion of the available cold gas sup- 
ply. The potential effect of such strong galaxy-galaxy inter- 
actions is of interest to studies of the environmental depen- 
dence of galaxy evolution as the merging rate (the product of 
relative velocity and interaction cross section) is predicted 
to be a strong function of environment: the cross section 
for disruptive encounters increases as the cluster velocity 
dispersi on approaches the in ternal velocity of the infalling 
galaxy (|Makino fc Hull 1 19971 '). Compelling evidence for en- 
hanced galaxy-galaxy interactions in rich cluster environ- 
ments is provided by Hubble Space Telescope (HST) ob- 
servations indicating a high merger fraction i n such environ- 
ments compared to field compa rison samples ([Dressier et al.l 
1 1994 1 van Dokkum et EtUl 19991 ). However, such observations 
must be contrasted with mid-infrared selected moderate 
starburst galaxies located in rich cluster environments whose 
optical morpholog i es resemble undistu rbed spiral galaxies 
jGeach et al.ll2009l : lOemler et a"Lll2009l ). 



A picture is therefore emerging whereby multiple mech- 
anisms (ICM stripping, merger induced star formation and 
tidally induced star formation) may participate in the phys- 
ical processing of infalling galaxies. Currently unanswered 
questions focus upon whether more than one physical pro- 
cess acts upon a typical galaxy falling into a dense envi- 
ronment (and which might be considered dominant) and 
whether the relative importance of each of the suggested 
physical processes changes as a function of the global proper- 
ties of the group/cluster environment into which the galaxy 
falls. In this paper we attempt to answer the related ques- 
tion of whether the effectiveness of galaxy processing can 
be determined as a function of global cluster environment. 
O ur approach is to c ompute the blue fraction as defined 
bv lButcher fc OemleJ lfl984 ) for a sample of X-ray selected 
galaxy clusters and to determine whether the blue fraction 
displays a significant trend versus X-ray temperature (here 
employed as a proxy for the global cluster mass). 



Throughout this paper, values of S7m.o ~ 0.3, flA,o ~ 0.7 
and Hq = 70kms~^Mpc~^ are adopted for the present 
epoch cosmological parameters describing the evolution of 
a model Friedmann-Robertson- Walker universe. All magni- 
tude information is quoted using AB zero point values. 
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□ XMM Clusters 
A CCCP Clusters 



Figure 1. Redshift and temperature distribution of the hot, mid 
and cool clusters samples as defined in the text. Open triangles 
represent CCCP clusters and open squares represent XMM-LSS 
clusters. 



scatter observed in scaling relations such as the X-ray L — T 
relation. In this sense, the CCCP clusters would all be de- 
tected by X-ray images of the same quality as the XMM- 
LSS survey images (irrespective of how the CCCP clus- 
ter were originally selected). However, XMM-LSS lacks the 
areal coverage to detect such massive clusters which display 
low sky surface densities. Therefore, the combined XMM- 
LSS/CCCP sample is a representative, (largely) unbiased 
sample of X-ray clusters with temperatures 1 < fcr(keV) < 
12. We note that the XMM LSS clusters populating the cool 
sample are moderately biased toward higher X-ray luminosi- 
ties than the average clu ster population at that temperature 
(see IPacaud et al.l 20071 Figure 8) . This implies a moderate 
bias toward gas rich or centrally condensed systems which 
we will recall at relevent points in the following discussion. 

The characteristic cluster scaling radii employed in 
this work are based upon rsoo. This is defined as the ra- 
dius at which the enclosed cluster mass density equals 
500 times the critical density of the universe at the 
cluster redshift JPacaud et al. 2007t ). Converting the rela- 
tion of iFi noguenov et al.| ( 200ll ') to a ACDM cosmology, 
IWillis et a l. (2005) define 



rsoo 



0.375 T°-^^h73{zy''Mpc 



(1) 



where T measured X-ray temperature in keV and /173 is the 
Hubble constant in units of 73 kms~^Mpc~^. 



2 THE GALAXY CLUSTER SAMPLE 

The data presented in this paper are drawn from two com- 
plementary samples of X-ray selected galaxy clusters. Clus- 
ters with X-ray temperatures T < 3 keV are drawn from 
the X-ray Multi-Mirror (XMM) Large Scale Structure (LSS) 
surve y. Clusters are selecte d from the 5 deg^ "Class 1" sam- 
ple of HacaudeFaD (|2007l ) according to X-ray temperature 
T < 3 keV and spectroscopic redshift 0.25 < z < 0.35 (Fig- 
ure [l}. Such relatively cool X-ray clusters are often referred 
to as galaxy groups in the literature. However, all X-ray se- 
lected systems considered in this paper are referred to as 
"clusters" for simplicity. The above criteria generate a sam- 
ple of 11 clusters which are referred to as "cool" in the fol- 
lowing analysis. 

Clusters with X-ray temperatures T > 5 keV are 
selected from the sample of Horner (2001) and form 
part of the Canadi an Cluster Comparison Project (CCCP; 
iBildfell et al]|2008l '). Although the final CCCP sample wiU 
contain some 50 X-ray selected galaxy clusters, we consider 
here the 27 systems with accompanying optical data ob- 
tained using the Canada France Hawaii Telescope (CFHT) 
MegaCam imager (see below). Due to the large temperature 
range covered by the CCCP sample (5 < T(keV) < 12) we 
further subdivide the CCCP sample into "mid" clusters dis- 
playing 5 < T(keV) < 8 and "hot" displaying T > 8 keV. 
We limit the CCCP sample to redshifts 0.15 < 2 < 0.4 in 
order to maximise the available sample size while reducing 
the redshift interval over which photometric quantities are 
fc-corrected to a common epoch. These considerations limit 
the size of the mid and hot samples to 18 and 5 clusters 
respectively. The X-ray properties of all clusters are shown 
in Table n 

The CCCP clusters all lie above the nominal X-ray fiux 
limit of the XMM-LSS survey and cover the full range of 



2.1 Optical photometry 

Optical photometry for the XMM-LSS and CCCP cluster 
samples is computed from CFHT Megacam images avail- 
able for all fields. XMM-LSS clusters lie within the CFHT 
Legacy Survey (CFHTLS) wide synoptic survej0 Wl area. 
The survey data consists of images taken in the CFHT ugriz 
filter set. Approximately 4 deg^ of the XMM-LSS footprint 
lies beyond the northern declination limit of the CFHTLS 
Wl, of which 3 deg'^ has been imaged as part of the XMM- 
LSS follow-up campaign using CFHT Megacam in the grz 
bands. Image exposure times in these additional fields are 
matched to the CFHTLS wide exposure times for each filter. 
The CCCP clusters considered in this sample were observed 
using CFHT Megacam as part of an optical follow-up of X- 
ray selected galaxy clusters (e.g. Bildfell et al. 2008). The 
processing of the CFHTLS wide data plus the northern ex- 
tension and of the CCCP Megacam data are described in 
Hoekstra et al. (2006) and Bildfell et al. (2008) respectively. 
Image data is available in the g and r bands. Table [5] de- 
scribes the main characteristics of each data set with Figure 
[5]showing r-band number counts for the CCCP and Wl data 
compared with the CFHT Deep survey. We adopt a magni- 
tude limit r = 23.5 in all subsequent analyses in order that 
photometric samples drawn from all three data sets can be 
reasonably considered to be complete. Given the overlap in 
the photometric filters employed by XMM-LSS and CCCP 
we consider the photometric properties of cluster galaxies in 
the combined sample as determined using g and r bands. 

Source extraction and photometry were performed us- 
ing SExtractor v2.5.0 (Berlin & Arnouts 1996). Image re- 
gions affected by saturated stars and detector artefacts were 

2 http : / /www.cfht . hawaii.edu /Science/CFHTLS / 
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Table 1. Properties of the cluster sample. Clusters are sorted 
with increasing temperature. Clusters possessing T(keV) < 3 are 
labelled "Cool", clusters possessing 5 < T(keV) < 8 are labelled 
"Mid" and those possessing T(kcV) > 8 are labelled "Hot" . 



Cluster 




R.A. 


Dec. 


Tx 


z 








(dcg.) 


(deg.) 


(kcV) 




(kpc) 


XLSSC 


13 


36.858 


-4.538 




0.31 


340 


XLSSC 


51 


36.498 


-2.826 


1 90-1-0.12 
^•^^-0.13 


0.28 


384 


XLSSC 


44 


36.141 


-4.234 


-1 97-1-0.09 
^■^'-0.1 


0.26 


399 


XLSSC 


08 


36.337 


-3.801 


1 30+"-^ 
l.ou_p 22 


0.30 


387 


XLSSC 


40 


35.523 


-4.546 


-0.25 


0.32 


402 


XLSSC 


23 


35.189 


-3.433 


1 57-1-0.29 
^■°'-0.16 


0.33 


497 


XLSSC 


22 


36.916 


-4.857 




0.29 


472 


XLSSC 


25 


36.353 


-4.679 


2 00+''''^2 


0.26 


533 


XLSSC 


18 


36.008 


-5.090 


2 n4+°-^^ 


0.32 


615 


XLSSC 


27 


37.014 


-4.851 


2 84+"-^ 


0.29 


653 


XLSSC 


10 


36.843 


-3.362 


2 40+" *^ 
'^■^"-0.36 


0.34 


574 


MS0440+02 


70.805 


2.166 


5 02+°-^l 


0.19 


957 


A1942 




219.600 


3.669 


5 12+°-''l 


0.22 


957 


A0223 




24.477 


-12.815 


5 28+°-''^ 

"^■^"-0.52 


0.21 


964 


A2259 




260.033 


27.668 


c- ,9-1-0.29 
"^■■^^-0.27 


0.16 


1007 


A1246 




170.9972 


21.482 


6 04+° *^ 


0.19 


1078 


A2537 




347.092 


-2.187 




0.30 


1039 


A0959 




154.433 


59.556 




0.29 


1022 


A0586 




113.072 


31.637 




0.17 


1127 


A0115 




13.980 


26.422 


p. AK-i-o.aa 

O-''3_0.31 


0.20 


1120 


A0611 




120.228 


36.065 


6-69+0.44 


0.29 


1100 


A0521 




73.510 


-10.244 


6.74l°i5 


0.25 


1106 


A2261 




260.609 


32.139 


6.88t°:ll 


0.22 


1153 


A2204 




248.192 


5.574 




0.15 


1247 


MS1008-12 


152.632 


-12.652 


7 -7-1-1.56 
' -1.21 


0.30 


1169 


CL1938+54 


294.555 


54.159 


y c:9-l-0.38 


0.26 


1200 


A0520 




73.554 


2.924 


7 01-1-0.74 


0.20 


1262 


A1758 




203.205 


50.538 


7 95+°-^* 

' ■^'-'-0.62 


0.28 


1186 


A2111 




234.914 


34.429 


8.02tO-?5 


0.23 


1267 


A0851 




145.746 


46.9945 


8-99ll:oi 


0.41 


1291 


A0697 




130.740 


36.3662 


14+0-6 
^•1^-0.54 


0.28 


1343 


A2104 




235.029 


-3.3017 


9 31+0-5 


0.16 


1438 


A1914 




216.512 


37.8244 


9 48+° *^ 


0.17 


1444 


A2163 




243.920 


-6.1438 


-19 -19+0.62 
-^^•-1^-0.57 


0.20 


1663 




Figure 2. Number counts as a function of 3" diameter r-band 
magnitudes in representative CCCP and CFHTLS Wide Mega- 
cam fields arc compared to number counts in representative 
CFHTLS Deep Megacam fields. The vertical dashed line indi- 
cates r = 23.5 and is the faintest magnitude employed in the 
following analysis. 



tude distribution for sources in each Megacam field (Fig- 
ure [Sj. The properties of instrument point spread function 
(PSF) that determines the half-hght radius of the steUar 
locus varies systematically over the Megacam field. Consid- 
ered over the entire field, this variation broadens the stellar 
locus and reduces the effectiveness of a single HLR cut at 
excluding stellar sources. The field itself is mapped onto a 
4x9 array of CCD detectors. We therefore determined the 
weighted mean HLR of the stellar locus in each Megacam 
CCD and rescaled the HLR values of all sources in each CCD 
such that the stellar values mapped on to a single locus de- 
fined in a reference CCD close to the optical centre of the 
detector. This operation reduced the effective dispersion of 
sources defining the stellar locus in each field (Figure |3]) and 
permitted the application of a single threshold to exclude 
stellar source^. 



Table 2. Characteristics of the optical data. 



Sample 


g-band 


r-band 




texp (s) 


te^p(s) seeing (") 


XMM-LSS 


2500 


2000 0.8 


CCCP 


1800 


4800 0.7 



excluded. Zero point information for sources detected in the 
CFHTLS Wl area plus northern extension was extrapolated 
from common sources detected in the Sloan Digital Sky Sur- 
vey equatorial patch which overlaps the southern edge of the 
Wl area. Source photometry is quoted in AB magnitudes 
measured within 3" diameter circular apertures. 

Star-galaxy separation was performed by considering 
the r-band half-light radius (HLR) versus r-band magni- 



3 COLOUR MAGNITUDE DIAGRAMS 

Colour magnitude diagrams for each of the hot, mid and 
cool cluster samples are displayed in Figures [S] [6] and [7] re- 
spectively. All sources lie within rsoo of the measured X-ray 
centre of each cluster. 

Computation of the blue fraction in each cluster first 
requires the red sequence relation defining the linear colour 
sequence followed by cluster early-type galaxies to be de- 
termined. This provides a reference value relative to which 

We note that a component of the PSF variation over the Mega- 
cam field is anisotropic. However, for the purpose of star-galaxy 
separation, application of an isotropic scaling factor is sufficient. 
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Figure 3. Location of the Stellar Locus in a typical CFHTLS Figure 4. Red edge diagram for XLSSC 22. See text for more 

Wide Megacam field, (a) Before PSF scaling is applied, (b) After details. 
PSF scaling is applied. See text for more details. 



the blue galaxy population in each cluster may be defined. 
To determine the location of the red sequence we first ap- 
plied a statistical background subtraction procedure to each 
cluster to highlight the population of c l uster members. We 
applied the method of IPimbblet et aP ()2002h whereby the 
field population and the cluster plus field population were 
each represented on a colour-magnitude grid. The probabil- 
ity that a galaxy occupying a particular grid position is field 
galaxy is then 

r,r 1 \ A X N{col,mag)Fieid 
P(col,mag)F,eid = -tt^. — , (2) 

Jy(cOl, magjciuster+Field 

were A is an areal scaling factor used to match the area of 
the field sample to the cluster area. The field colur magni- 
tude distribution is computed employing all galaxy sources 
at radii > 8 x rsoo from the X-ray cluster centre. For each 
galaxy within rsoo of the cluster centre, membership was 
determined by comparing a random number in the inter- 
val [0,1] to the field probability value. This procedure was 
repeated 100 times for each cluster. 

For each realisation, the red sequence location was com- 
puted employing a weighted, linear least squares fit and the 
mean slope and zero-point was calculated from the distri- 
bution of 100 values. To check and, if necessary, refine the 
location of the mean red sequence for each cluster we consid- 
ered the "red edge" diagram for each system (an example of 
which is displayed in Figure The red edge diagram shows 
the number of galaxies located at a colour offset A{g~r) at 
fixed r-magnitude from the red sequence in a particular clus- 
ter. The zero-point of the red sequence was then adjusted 
to set the peak in the red edge distribution at zero colour 
offset. 

This procedure worked well for the mid and hot cluster 
samples. However, following statistical background subtrac- 
tion, the cool clusters typically possess insufficient cluster 



members to generate a reliable red sequence fit. In this case 
the red sequence relation computed for the stacked mid clus- 
ter sample described in Section [5] was applied to the cool 
clusters. The red sequence zero-point was adjusted using 
the red edge diagram procedure to provide an improved fit 
to individual clusters. 

Figure |8] indicates the g — r colour of the red sequence 
measured at fixed absolute magnitude (My = —21) versus 
redshift for the mid and hot clusters. The colour uncertainty 
is computed as the standard deviation of the red sequence 
zero-points as measured for the 100 realisations of the cluster 
subtraction procedure. The fc-correction employed to con- 
vert absolute V-magnitude to apparent r-magnitude at the 
redshift of each cluster assumes a early-type galaxy spectral 
energy distribution (SED; Kinney et al. 1996). We note that 
the exact choice of template is not a significant factor when 
computing the apparent r-magnitude reference location on 
the red sequence as the fc-correction is dominated by the 
bandwidth term at 2 ~ 0.3. Furthermore, due to the small 
slope of the fitted red sequence relation in each case ('^ 0.05) 
small systematic magnitude errors result in negligible colour 
uncertainties. 

The best fitting SED describing the observed colour evo- 
lution of the mid and hot cluster red sequences is computed 
employing a linear interpolation between an early-type (Ell) 
and early- type spiral (Sab) SED, e.g. 

SED(A) = (1 -s) EU(A) -f a; Sab(A), (3) 

where the parameter x is computed using a minimisa- 
tion procedure. The best fitting SED model was found to 
be approximately 70% Elliptical and 30% Sab. We employ 
this template as a reference point from which to determine 
accurate fc-corrections and blue galaxy colour thresholds for 
all clusters (cool, mid and hot) in the following discussion. 
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Figure 5. Hot sample clusters. All sources within rsoo of the cluster centre are plotted. The red 
line marks the location of the red sequence and the blue line marks the Butcher & Oemler (1984) 
cut as described in the text. The vertical dashed line indicates the r magnitude corresponding to 
My = — 20 at the cluster redshift. 



4 BLUE FRACTIONS 

The blue fraction of each cluster was computed foUowing 
the definition of Butcher & Oemler (1984). Following their 
approach, all galaxies displaying an absolute magnitude 
My ^ — 20 are considered and blue galaxies are defined as 
those displaying a rest frame colour offset A(_B — V) = —0.2 
measured relative to the red sequence. We add the further 
criterion that galaxies within rsoo of the cluster X-ray cen- 
tre are considered part of the "total" , i.e. cluster plus field, 
populations, and galaxies at clustercentric radii > 8 x rsoo 
and within the same Megacam field are considered as the 
"field" population. The blue fraction within each cluster is 
then computed as 

J _ Nsiue.Total — A NbIuc, Field ^^-^ 



where NBiue,Totai is the number of blue galaxies in the clus- 
ter plus field, Nsiue, Field IS the number of blue field galax- 
ies, Nrotai is the total number of galaxies in the cluster plus 
field and NFieid is the total number of field galaxies. The 
symbol A denotes an areal scaling factor to correct the field 
population area to that of the cluster area. 

The rest frame Butcher & Oemler (1984) magnitude 
and colour criteria were expressed as observed frame r- 
magnitude and A((/ — r) colour offsets at the redshift of 
each cluster by considering the galaxy SED implied by the 
cluster red sequence colour (Section [31). The colour offset 
was computed by generating a second SED model according 
to Equation [3] displaying A(_B ~V) = —0.2 compared to the 
reference SED describing the cluster red sequence relation. 
This "blue-cut" SED model was found to be 75% Sab and 
25% Sbc. Employing this SED model A((? — r) correspond- 
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Figure 6. Mid sample clusters. All sources within rsoo of the cluster centre are plotted. The red 
line marks the location of the red sequence and the blue line marks the Butcher & Oemler (1984) 
cut as described in the text. The vertical dashed line indicates the r magnitude corresponding to 
My = — 20 at the cluster redshift. 



ing to A{B ~ V) — —0.2 was calculated for each cluster at 
the appropriate redshift. 

The blue fraction error is estimated by computing a dis- 
tribution of blue fraction values for each cluster. Each blue 
fraction value is computed as follows: circular apertures of 
radius rsoo are placed at random locations within the Mega- 
cam field of each cluster. The galaxy population within each 
aperture is employed as the background value for Equation 
|4] The error on the blue fraction for each cluster is then 
estimated as the interval about the median blue fraction 
value containing 67% of the distribution. Errors computed 
using this method are typically 1.7 times larger than those 
computed assuming Poissonian uncertainties alone. 

Figure [9] displays the blue fraction as a function of red- 
shift for all clusters and shows an apparent trend to observe 
increasing blue fraction versus redshift. However, splitting 
the sample by X-ray temperature reveals that this trend 



may instead arise from the varying global environment of 
each cluster modulo the slightly different redshift interval 
covered by each of the cool, mid and hot samples. 

Figure [TD] displays the blue fraction as a function of 
cluster X-ray temperature for the cool, mid and hot samples. 
In addition, blue fraction values and Poisson uncertainties 
are listed in Table O The data indicate that typical blue 
fraction in each cluster and the dispersion in blue fraction 
values within a given temperature sub-sample increase as 
the temperature of the X-ray cluster decreases. 

Cluster s of s imilar properties have been studied by 
IWake et al.l (|2005l ) and we also compare their results to ours 
as a function of cluster X-ray temperature in Figure 1101 
One immediately notes that the blue fraction values for the 
current sample are consistently larger than the sample of 
Wake et al. though nominally covering a similar range of 
X-ray temperature. Blue fractions in the Wake et al. sample 
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Figure 7. Cool sample clusters. All sources within rsoo of the cluster centre are plotted. The red 
line marks the location of the red sequence and the blue line marks the Butcher & Oemler (1984) 
cut as described in the text. The vertical dashed line indicates the r magnitude corresponding to 
Mv = —20 at the cluster redshift. The triangles represent spectroscopically confirmed members. 



are computed within an aperture of radius one-third of the 
virial radius. This corresponds to a radius approximately 
equal to 0.7 rsoo and we have corrected the Wake et al. 
values to the aperture used in the current paper using the 
blue fraction versus radius curve appropriate to each clus- 
ter temperature (shown in Figure I14[) . This correction is 
necessarily approximate and clearly offsets between the two 
samples remain. The remaining differences arise from the 
different methods used to determine blue fraction values in 
each sample. Specifically, Wake et al. assume that the SED 
representing the cluster colour magnitude relation is a pure 
elliptical model rather than the composite elliptical plus Sab 
model adopted in the current paper. The use of different 
spectral models results in different fe-corrections versus red- 
shift with the consequence that the observed frame colour 
offset employed to define to the blue cut in each cluster is 
typically larger in the Wake et al. sample than that employed 



here. The application of a bluer cut relative the OMR in a 
given cluster results in a smaller blue fraction. We therefore 
conclude that the remaining offset in blue fraction versus 
temperature between the Wake et al. sample and that pre- 
sented in the current paper results from such methodology 

differences. 

As noted bv lMargoniner et al.l l|200ll ) and [Hansen et ahl 

(|2009l ). the observed trends in cluster blue fractions can be 
explained as a function of both redshift and cluster mass. In 
order to determine the relative influence of varying redshift 
and cluster temperature (here used as a proxy for cluster 
mass) upon the blue fractions we flt the data for the three 
samples with a simple function of the form 

fB{z,T)^Po + P.z + l3TT-^ (5) 

where Po, 13 z and /3t are constants to be determined employ- 
ing a minimum algorithm and T is expressed in keV. The 



An environmental Butcher- Oemler effect in intermediate redshift X-ray clusters 9 





5 10 

X-ray temperature (keV) 



Figure 8. The red crosses indicate the g — r colour of the fitted 
red sequence relation measured at My = — 20 for the mid and 
hot samples. The dashed line marks the 100% elliptical model for 
reference. Blue crosses indicate the location of the red sequence 
in the cool clusters determined using the mid cluster red sequence 
relation (sec text for more details). 




Figure 9. Cluster blue fractions as a function of redshift. 



best fitting values are = 0.04 ±0.02, = 0.67 ±0.08 and 
Pt = 0.42 ±0.07. Confidence intervals of the fitted values of 

and Pt are shown in Figure [11] The minimum value of 
obtained using this procedure is approximately 3 per degree 
of freedom and, if the blue fraction errors are accurate, may 
indicate the presence of intrinsic scatter in the distribution 
of blue fraction values. The current sample of XMM-LSS 



Figure 10. Cluster blue fraction as a function of X-Ray temper- 
ature. 



and CCCP clusters is insufficiently large and does not cover 
a large enough interval in either redshift or temperature to 
constrain both the exponents of the redshift and tempera- 
ture dependence of Equation[5]in addition to the coefficients. 
We can rule out no redshift evolution but the temperature 
dependence is relatively unconstrained (and is complicated 
by the degeneracy with /3o as the power law exponent of T 
tends to zero). 

An alternative approach is to investigate the trends in 
the data employing a non-parametric approach. In this case 
we apply a partial Spearman rank analysis to the data in or- 
der to examine the extent of any correlation between /s and 
either temperature or redshift while controlling for the varia- 
tion of the second variable. The partial correlation coefficient 
describing blue fraction and redshift while controlling for 
temperature is r(/s,2,T) = 0.39 while the corresponding 
coefficient describing blue fraction versus temperature while 
controlling for redshift is r{fB,T,z) = —0.64. The proba- 
bility of obtaining r — —0.64 from the sample by chance is 
approximately 0.2% whereas the corresponding probability 
for the value r = 0.39 is somewhat greater at approximately 
1.8%. Overall, both the parametric and non-parametric ap- 
proach confirm the evolution in the cluster blue fraction with 
both redshift and cluster temperature. The significance of 
the greater blue fraction versus decreasing temperature can 
also be seen in Figure [13] where the average blue fraction for 
the cool sample is significantly higher than either the blue 
fraction for the mid and hot samples. 

One may further investigate the dependence of blue 
fraction on cluster temperature by correcting /s values for 
individual clusters to a common epoch a,t z — 0.3 using 
Equation [5] and the best-fitting coefficients. The corrected 
fs values are displayed in Figure [12] and confirm the trend 
to observe greater blue fraction values in cooler clusters. We 
consider further the implications of this trend in Section [6] 
We next consider the effects of the various choices made 
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Figure 12. Cluster blue fraction as a function of X-Ray tem- 
perature. The black points indicate the result of correcting the 
original blue fraction values to a common epoch ai z = 0.3. 

during the blue fraction analysis on the overall robustness 
of the blue fraction trends versus redshift and temperature. 

4.1 Testing the blue fraction computation 
assumptions 

The cluster blue fraction may be computed employing one 
of a number of crtieria to segregate the red and blue clus- 
ter galaxy populations. We apply the original definition of 
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Figure 13. A comparison of blue fraction computation methods 
versus temperature. The circular points indicate the average blue 
fraction in each of the three temperature sub-samples computed 
using the B084 method. The triangles indicate the average blue 
fraction per temperature sub-sample computed using the defini- 
tion f'g = 1 ~ fii (see text for details). The triangles have been 
offset in temperature from the circles for clarity 



iButcher fc Oemleil dlOsil ) in this paper but note the critis- 
cisms of lAndreonI (|2005l ) who argue that adopting a fixed 
rest frame colour offset relative to the red sequence is prone 
to give misleading blue fraction trends as galaxies of dif- 
ferent SEDs drift in and out of the adopted colour interval 
as a funtion of redshift. To assess whether this is an im- 
portant consideration in our analysis we compare the aver- 
age blue fraction per cluster sub-sample to a modified blue 
fraction determined by fitting the red fraction associated 
with the stacked CMD of clu sters by temperat ure. We fol- 
low the method described by iLoh et ahl l|2008l ) and model 
the red wing of the red sequence of the CMD of all clus- 
ters stacked by temperature once fc-corrected to a common 
redshift z — 0.3 (we defer the description of the stacking pro- 
cedure to Section [Sjl. The red wing is modelled as a double 
Gaussian function and the best-fitting function is reflected 
about the location of the red sequence to determine the num- 
ber of red galaxies in each distribution. We then compute 
the red fraction in an analogous manner to the blue fraction 
and define a modified blue fraction as f'g — 1 — fa- We com- 
pare the original and modified blue fraction values for each 
temperature sub-sample in Figure [13] and note that each 
approach reveals the same trend (and that the points are 
identical within the errors). We are therefore satisfied that 
the original B084 definition produces reliable results and we 
avoid the requirement to fit double Gaussian features to in- 
dividual cluster colour distributions which, in the case of the 
cool clusters, typically result in poorly constrained models. 

We next consider the variation of the computed blue 
fraction with assumed cluster radius to a) investigate 
whether the choice of radius introduces a bias into the trend 
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Table 3. Blue Fractions for all clusters employing a limiting mag- 



nitudc My = 


—20 within rsoo 


of the cluster centre. 


Cluster 


/s 


Number of blue galaxies 


XLSSC 13 


0.69 ± 0.20 


11 


XLSSC 51 


0.55 ± 0.15 


11 


XLSSC 44 


0.61 ± 0.14 


13 


XLSSC 08 


0.50 ±0.23 


14 


XLSCC 40 


0.40 ±0.17 


5 


XLSSC 23 


0.35 ± 0.05 


16 


XLSSC 22 


0.53 ±0.15 


10 


XLSSC 25 


0.56 ±0.10 


17 


XLSSC 18 


0.44 ± 0.09 


14 


XLSSC 27 


0.59 ±0.12 


12 


XLSSC 10 


0.53 ±0.12 


20 


MS0440+02 


0.19 ± 0.07 


16 


A1942 


0.31 ± 0.U4 


An 
49 


AU223 


ri oo 1 n r\ A 

0.23 ± 0.U4 


30 


A2259 


0.40 ± 0.U5 


50 


A 1246 


nod 1 n nti 

0.24 ± 0.U2 


30 


A2537 


0.31 ± O.Ub 


54 


A0959 


0.21 ± 0.05 


25 


A0586 


0.22 ± 0.04 


39 


A0115 


0.25 ± 0.04 


38 


A0611 


0.36 ±0.03 


65 


A0521 


0.24 ±0.05 


55 


A2261 


0.38 ±0.05 


80 


A2204 


0.31 ± 0.03 


35 


MS1008-12 


0.27 ±0.05 


58 


CL1938-I-54 


0.38 ±0.03 


95 


A0520 


0.21 ± 0.02 


43 


A1758 


0.23 ±0.03 


52 


A2111 


0.15 ± 0.02 


26 


A0851 


0.38 ± 0.02 


97 


A0697 


0.24 ±0.02 


54 


A2104 


0.14 ±0.02 


24 


A1914 


0.18 ± 0.03 


35 


A2163 


0.23 ± 0.03 


69 



of observed blue fraction versus temperature, b) determine 
the effect of uncertainty in the estimated rsoo value for each 
cluster on fs and c) relate blue fractions computed for the 
current sample to comparable studies in the literature. Fig- 
ure [13] displays the average blue fraction per temperature 
sub-sample computed within circular apertures of radius ex- 
pressed as a fraction of rsoo. This variation of the blue frac- 
tion with characteristic radius has been noted in previous 
studies and th e trend displayed in Figure | 14 | are consistent 
wit h those of iFairlev et all l|200ll) . lEUingson erahl (120011 ) 
and I Wake et al.l ( 2003)- Though we discuss this trend further 
in Section [6] we note here that the increasing trend of blue 
fraction with scaled aperture radius is nominally consistent 
with a blue infalling field population that is processed by 
ram pressure stripping upon falling into each cluster. How- 
ever, the details of this conclusion will be discussed in Sec- 
tion[6] The average blue fraction at rsoo in each temperature 
sub-sample reflects the trend observed in Figure [TO] What 
can be noted at this point is that the exact choice of radius 
wihtin which blue fractions are computed does not affect the 
conclusion that cooler clusters display greater blue fractions 




Figure 14. Radial variation of the cluster blue fraction 

than hotter clusters (within the range 0.5 rsoo < r < 2 rsoo 
investigated herej3. 

Uncertainty in the temperature computed for each clus- 
ter propagates to an error in the value of rsoo computed for 
each cluster. The typical temperature error for the XMM- 
LSS clusters contributing to the cool sample are of the order 
of 10-20% (Willis et al. 2005). The gradient of the trend of 
/s versus r/rsoo estimate at rsoo is approximately 0.15. This 
indicates that uncertainty in cluster X-ray temperature will 
contribute a fractional fs error of order a few percent, i.e. 
small compared to the Poissonian error in fs computed for 
an individual cluster. 

The mid temperature CCCP cluster MS 1008 i s coni - 
mon to both this work and that of lKodama fc Bower! (1200 ll ) 
and provides a useful check upon our results. The value 
of the blue fraction obtained by these authors was fn = 
0.161±0.036 and was calculated within R30 (the radius con- 
taining 30% of the total number of galaxies in the cluster) 
corresponding to an angular radius of 2.82'. We obtain a 
value /s — 0.27±0.04 within rsoo corresp onding to an angu- 
lar rad ius of 4.44'. The radial cut used bv lKodama fc Boweij 
('2001') corresponds to approximately 0.63 rsoo and when the 
blue fraction was recalculated within this fraction of rsoo, 
a value of /s = 0.18 ± 0.04 was obtained, improving the 
agreement with lKodama fc Boweij (|200ll ). We note that the 
errors quoted here are purely Poissonian and do not account 
for the effect that minor differences in the methodology used 



^ Furthermore, we note that the value of rsoo in each case is 
computed employing the observed cluster temperature and an as- 
sumed mass-temperature relation. X-ray spectral measurements 
from which the cluster temperature is computed are typically ex- 
tracted within an on-sky aperture of size of the order of rsoo in 
each cluster. Therefore, employing computed rsoo values to define 
aperture sizes up to several times the virial radius would invlove a 
considerable extrapolation beyond the scale on which the X-rays 
are measured. 
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Figure 15. Variation of the cluster blue fraction for each cluster 
sample with faint magnitude cut. 

in each study (e.g. fe-correction, computation of the red se- 
quence relation) have upon the computed value of /s. 

We next consider the effect of extending the faint abso- 
lute magnitude cut applied when computing the blue frac- 
tion. Figure [15] displays the average blue fraction per tem- 
perature sub-sample computed within a circular aperture 
of rsoo as the applied faint magnitude cut is varied from 
Mv = —22 to -18. The main feature of the diagram is that 
the excess fraction of blue galaxies observed in cool clusters 
compared to the mid and hot samples continues to increase 
as the sample magnitude limit is extended to fainter magni- 
tudes. 

When considered together, the variation of blue frac- 
tion with both radius and magnitude indicate that the ex- 
cess blue galaxy fraction observed in cool clusters is domi- 
nated by the contrubtion of faint galaxies at large radius. In 
the following section we attempt to identify this population 
directly on the colour-magnitude diagram for each temper- 
ature sub-sample. 



5 STACKED COLOUR MAGNITUDE 
DIAGRAMS 

The previous analyses have identified a relative excess in 
the fraction of blue galaxies in cool X-ray clusters com- 
pared to the mid and hot samples. While the blue frac- 
tion is a useful quantity to investigate, it may prove infor- 
mative to consider the average CMD for each temperature 
sample for which blue fractions are computed. The statis- 
tical background subtraction method described in Section 
[3] was applied 100 times to each cluster and the average 
CMD was computed by stacking the individual CMDs on 
a binned colour magnitude plane in intervals of 1 mag. in 
magnitude and 0.5 mag. in colour. The average background 
subtracted CMD for each cluster was then transformed to 



a common redshift z = 0.3. The transformation described 
the effects of distance dimming and the fe-correction. The fc- 
correction for each colour pixel on the CMD plane was com- 
puted using the best-fitting interpolated spectral template 
required to reproduce the required observed colour value at 
the cluster redshift. Individual cluster CMDs within each 
temperature sub-sample were then stacked on the z = 0.3 
CMD plane. Each cluster CMD was assigned equal weight in 
the stacking process by normalising all transformed CMDs 
to a total contribution of unity summed withing the region 
21 < r < 22 and 1 < g ~ r < 1.5. The stacked colour mag- 
nitude diagrams for each sample are shown in Figure 1161 
In addition, the fourth panel in Figure [16] indicates the re- 
sult of subtracting the mid cluster stacked CMD from the 
corresponding cool cluster sample to highlight the location 
on the CMD of the excess fraction of blue galaxies found 
in the cool sample. As indicated in the plots, the relative 
excess appears uniform across a range of colours bluer than 
the red sequence. There is some evidence for an excess of 
faint (r > 23) blue (g — r ~ 0.2) galaxies in the cool sample. 
However, the CFHTLS optical data used to construct the 
cool CMD begins to be incomplete at these magnitudes and 
the purely visual impression may be flawed. 

An alternative approach is to consider the colour distri- 
bution of all galaxies brighter than some limit. The colour 
distributions derived from the unweighted sum of all clus- 
ters in each temperature sample (corrected to z = 0.3) are 
shown in the upper three panels of Figure 1171 The limiting 
magnitude is r — 21.5, approximately equal to Mv = —20 at 
this redshift. In order to investigate the distribution of blue 
galaxies in each sample the red s equence wa s mod elled and 
removed in each case. Following iLoh et al.l (120081 ) the red 
wing of the red sequence was fitted using a double Gaussian 
model with the mean of each Gaussian set to the location 
of the red sequence and the best fitting full-width at half- 
maximum (FWHM) and normalisation of each profile deter- 
mined using a minimum algorithm. The resulting model 
is then reflected about the location of the red sequence and 
subtracted from the corresponding colour distribution. 

The mean g — r colour of all galaxies bluer than the 
location of red sequence for each sample was calculated and 
found to be 0.70, 0.675 and 0.675 for the cool, mid and hot 
samples respectively. We hesitate to place too much empha- 
sis on the trend to observe redder mean colours in the blue 
cloud for lower temperature X-ray systems. This is princi- 
pally due to potential limitations such as a) the obvious sub- 
traction artefact in the mid sample red sequence subtracted 
distribution and b) the broad nature of such peaks and the 
low number of galxies involved (in the cool systems). 

The red sequence subtracted distribution for the cool 
sample does appear to show a population of extremely blue 
galaxies not seen in the corresponding mid and hot cluster 
samples. To quantify this additional population we define 
the extremely blue fraction to be the fraction of galaxies 
bluer than the colour of an Sc galaxy {g — r = 0.257 at 
z — 0.3) for all galaxies down to r = 21.5 for each sample. 
These extreme blue fractions were found to 0.028 ± 0.004, 
0.0006 ± 0.0001 and 0.001 ± 0.0002 for the cool, mid and 
hot samples respectively. It was noted in Figure [15] that the 
fraction of blue galaxies in the cool cluster sample increases 
with the increasing faint magnitude limit used to define the 
blue fraction. With this in mind we computed the extremely 
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Figure 16. Stacked CMDs for the cool, mid, iiot and (cool-mid) cluster samples (see text for details). All data are 
fe-corrected to z = 0.3. In each panel the solid line marks the location of the red sequence and the dashed line marks 
the location of the corresponding Butcher-Oemler blue cut. 



blue fraction for the cool sample using galaxies in the mag- 
nitude range 21.5 < r < 23 as shown in Figure 1181 The 
computed extreme blue fraction is equal to 0.052 ± 0.004, 
i.e. the population of extremely blue galaxies increases with 
increasing magnitude. 



6 DISCUSSION 

We have presented an analysis of the fraction of bhie galax- 
ies in a large sample of X-ray clusters spanning a wide range 
of X-ray temperature with uniform optical photometry. We 
have noted that the computed blue fraction values display a 
trend to decrease with both decreasing redshift and increas- 
ing temperature. 

The trend of decreasing cluster blue fraction versus 
redshift has been observed previously and has been in- 
terpreted as a global effect driven by the decreasing rate 
of field galaxy infal l onto the cluster environment (e.g. 
lEUingson et al.ll200ll ). Coupled with the increasing cluster 
mass with decreasing redshift, the specific infall rate per 
unit cluster mass decreases at an even faster rate. Infall 
effects are further compounded with the sharply declining 
globa l star formation rate as a function of decreasing red- 
shift (|Hopk ms fc Beacomllioog ) which may be viewed as a 
decrease in the available gas supply of infalling galaxies. 
Looking beyond the redshift dependent component of the 
Butcher-Oemler effect, the decreasing blue fraction with in- 



creasing X-ray temperature identified in this paper is con- 
sistent with an environmental component to the Butcher- 
Oemler effect. 

The trend of decreasing blue fraction versus tempera- 
ture is most apparent whan comparing the across XMM-LSS 
and CCCP samples rather than within either sample. It re- 
mains possible that this difference is due in part to some 
subtle selection effect arising during the construction of the 
two samples rather than as a result of some physical pro- 
cess linked to the global environment. However, it is not 
clear what form such a bias would take. At this point we re- 
call that the cool clusters may display a mild bias to higher 
gas densities or more centrally concentrated systems com- 
pared to the average population. One would naively expect 
such systems to be more effective at processing infalling field 
galaxies compared to the average (i.e. lower X-ray luminos- 
ity) cluster over the same temperature interval - effectively 
introducing a mild bias toward lower blue fractions in this 
temperature sub-sample. One could therefore argue that ac- 
cess to a more complete census of the low-temperature clus- 
ter population would only amplify the trend of blue fraction 
versus temperature. It has also been noted that optically se- 
lected clusters display different optical properties compared 
to X-ray selected samples (see Haines et al. 2009 for a dis- 
cussion). However, X-ray samples of bright clusters are typi- 
cally complete in terms of the optical properties they sample 
(mainly because bright X-ray clusters are rare and thus it is 
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Figure 17. The top panels show stacked histograms for the cool, mid, and hot samples before red sequence sub- 
traction. All cluster data are background subtracted and fc-corrected to 2 = 0.3. The absolute numbers of galaxies 
in each bin have been re-scaled purely for visualisation purposes. The vertical line indicates the Butcher-Oemler 
blue cut location at z = 0.3. The bottom panels show the stacked histograms after red sequence subtraction (see 
text for details). The vertical lines incidate the observed frame colours of the Elliptical, Sa and Starburst 1 (SBl) 
models of Kinney et al. (1996). 



possible to select complete samples). Moreover, clusters se- 
lected by single band optical observations (i.e. purely based 
upon the projected overdensity of bright galaxies) are more 
likely to be biased toward high blue fraction values due to 
the high mass-to-light ratios of blue star forming galaxies. 
We therefore note that any possible bias between the XMM- 
LSS and CCCP samples would be naively be expected to 
generate the opposite trend to that observed in this paper 
and we remain confident that the trend of blue fraction ver- 
sus temperature we have identified is not a result of sample 
selection effects. 

A currently favoured explanation for the observed trend 
of blue fraction versus temperature is that infalling field 
galaxies are processed physically as they interact with the 
cluster environment. This interaction may take the form of 



ram pressure stripping whose effectiveness is a relatively 
simple function of the mass scale represented by the group or 
cluster into which the galaxies are falling. An alternative ex- 
planation is that infalling galaxies are processed via galaxy- 
galaxy interactions and therefore respond more readily to 
a combination of local rather than global velocity disper- 
sion and galaxy density. Whatever the cause, it is clear that 
the process by which blue galaxies are processed to appear 
red is more complete in hotter (more massive) environments 
compared to cooler (less massive) environments. 

Disentangling the above two effects is difficult as, to first 
order, the effects upon star formation in the infalling galaxy 
depend on environment in the same manner: ICM stripping 
and subsequent SF suppression are expected to be weakest 
(though not absent) in group environments - galaxy popu- 
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Figure 18. Red sequence subtracted histogram for thie cool sam- 
ple. Eacli iiistogram is computed ai z = 0.3. The upper panel dis- 
plays background correctd sources with r < 21. The lower panel 
displays background corrected sources with 21.5 < r < 23. The 
vertical lines incidate the observed frame colours of the Elliptical, 
Sa and Starburst 1 (SBl) models of Kinney et al. (1996). 



lations in these environments will display a larger fraction of 
blue galaxies relative to more massive clusters - while SF en- 
hancements associated with galaxy-galaxy interactions will 
be greater on group scales. Each effect would naively gener- 
ate the same observed trend of blue fraction versus temper- 
ature. 

The blue fraction in each temperature sub-sample dis- 
plays a very similar trend versus scaled radius and is consis- 
tent with an infalling field populatio n albeit the cool sam - 
ple is offset to higher blue fraction (|Ellingson et al.ll200l[ ). 
However, the cool sample displays a marked increase in blue 
fraction with decreasing magnitude. The increase in the cool 
sample blue fraction with increasing magnitude results from 
both a moderate increase in the number of blue galaxies but 
also in a relative deficit in the numbers of red galaxies com- 
pared to those in the mid and hot samples. We will address 
the red sequence luminosity function of the cool, mid and 
hot samples in a subsequent paper. However, we note that 
the relative deficit of red galaxies in cool clusters compared 
to the mid and hot samples is consistent with the scenario 
where the processing of blue galaxies to red galaxies during 
infall is relatively incomplete. 

Having subtracted the red sequence contribution from 
each colour histogram the blue galaxy colour distribution for 
each temperature sub-sample appears very similar. We do 



note an excess of extremely blue galaxies in the cool sample 
compared to the mid and hot samples. The excess is small 
(an extremely blue fraction of 3% compared to < 1%) yet 
significant. In addition the fraction of extremely blue galax- 
ies increases in the cool sample from ~ 3% to ~ 5% for 
galaxies displaying r < 21.5 and 21.5 < r < 23 respectively. 
One may tentatively associate these galaxies as actively 
star-bursting gal axies potentially driven by galaxy-galaxy 
intera ctions (e.g. iDe Propris et al.l [20031 ) . lOe Propris et al.l 
comment that the absence of a strong Butcher- 
Oemler effect in near-infrared selected cluster galaxy popula- 
tions compared to optically selected populations may be due 
to an increasing contribution from faint blue dwarfs whose 
optical brightness is boosted by recent star formation. The 
extremely blue population identified in the cool sample is 
nominally consistent with this explanation yet clearly forms 
only a small component of the overall blue cluster galaxy 
population. The small relative contribution of such poten- 
tially star bursting galaxies to the cluster population may 
be understood in terms of the short timescale over which a 
recent star burst will affect the integrated colour of an es- 
ta blished stellar popul ation. This issue has been addressed 
by iBarger et al.l (Il996l ) who simulate the colour evolution 
with time of an underlying passive (elliptical-type) and con- 
tinuous star formation (spiral-type) stellar population expe- 
riencing a short (~ 100 Myr), secondary (~ 10% by mass) 
burst of star formation associated with an interaction result- 
ing from the cluster environment. The shift towards bluer 
colours as a result of a moderate burst of star formation lasts 
typically only as long as the burst duration before the lumi- 
nosity weighted colour reverts to that of the underlying pre- 
burst stellar population. The effect is more pronounced when 
one applies the condition that all star formation ceases af- 
ter the secondary burst, in which case the continuously star 
forming stellar population (spiral-type) reddens to approx- 
imately the same colour as the passive population within 
~ 1 Gyr of the burst. 

Our analysis has identified a clear environmental depen- 
dence in the blue fraction of galaxies. However, the cause of 
this environmental Butcher-Oemler effect - whether it be 
ram pressure stripping or galaxy-galaxy interactions - can- 
not be discerned unambiguously on the basis of the colour 
distribution of the blue galaxy population in each sample. 
Clearly further information is required, most notably the 
combination of colour information with morphological infor- 
mation. One would naively expect the above two environ- 
mental processes to affect galaxy morphology in markedly 
different ways, e.g. ram pressure effects hould not disrupt 
galaxy disks whereas a strong disruptive effect is expected 
from galaxy-galaxy interactions. We will present this inves- 
tigation in a future paper. 
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